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ABSTRACT 


The earth's desert boundary layer is modeled in a 
wind tunnel. Velocity profile measurements are deter- 
mined on and around a barchan dune model inserted in the 
roughness layer on the tunnel floor. A theoretical in- 
vestigation is made into *che factors influencing the rate 
of sand flow around the dune. Flow visualization tech- 
niques are employed in the mapping of strec^mlines of flow 
on the dune's surface. Maps of erosion and deposition of 
sand are constructed for the barchan model, utilizing 
both flow visualization techniques and friction velocities 
calculated from the measured velocity profiles. The 
sediment budget found experimentally for the model is 
compared to predicted and observed results reported by 
others . The comparison shows fairly good agreement 
between the experimentally determined and predicted sedi- 
.ment budgets; this provides encouragement for the utili- 
zation of wind tunnel simulation and friction velocity 
approximation in determining erosion and deposition rates 
due to wind action. 
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Chapter 1 


INTRODUCTION AND BACKGROUND 

1.0 Introduction and Background 

Photographs taken from the Mariner 9 spacecraft of 
the surface of Mars reveal the presence of eolian (wind 
blown) landforms closely resembling the desert dunes 
found on Earth. In order to study the characteristics 
of the atmospheric winds and surface features of Mars, it 
was determined that more information v;as needed concern- 
ing the effects of the earth's desert boundary layer 
winds on the formation and transport of sand dunes. 
Measurements of the velocity profiles surrounding actual 
barchan dunes were made at a location near the Salton 
Sea in California [1] . The barchan dune, a crescent- 
shaped sand dune, was selected as an appropriate land- 
form for study because it appears in the Mariner 9 photo- 
graphs- It occurs as an isolated phenomenon and as the 
result of unidirectional winds convenient to simulate 
in a wind tunnel. The barchan is also self-preserving 
in shape and size over time. 

A wind tunnel study was performed, wherein the earth 

desert boundary layer was simulated, and scale models of 

the barchan dunes observed in the field were inserted. 

One of the mo'N^ls, selected on the basis that it was 

symmetrical and small enough to rotate through angles 

1 
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of -30® to 30® in the tunnel, was subjected to greater 
investigation than the others. Measurements of velocity 
profiles around the model yielded information about the 
local friction velocities on the dune, from which ex- 
pressions could be derived for sand transport on and 
around the dune. In addition, qualitative information 
about the wind flow at various heights above the surface 
of the dune was obtained from the profiles. Streamline 
maps determined from hotwire measurements, flow visuali- 
zation methods and measurements made in the field aided 
in predicting the sand flow. 

The "dust bowl" erosion of the Midwest in the 'thirties 
stimulated a number of agriculture-related studies of wind 
erosion, focusing on the characteristics of particle salta- 
tion and prevention of erosion. Saltation is a method of 
sand transport characterized by a jumping movement of 
sand particles. Semi-empirical formulae were developed 
to give the amount of sand transported over a flat surface 
as a function of grain size and friction velocity. 

Theoretical studies made of the forces on a sand 
grain in saltation were helpful in determining the factors 
to be considered in calculating the sediment budget for 
the barchan dune. Some controversy arises in determining 
the forces responsible for the initiation of saltation; 
in particular, it is not clear what role lift plays in 
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Experimentally determined maps of erosion and deposi- 
tion of sand on the barchan dune, constructed from stream- 
line maps and a friction velocity contour map, are com- 
pared to a predicted map of the sediment budget (for dune 
equilibrium) and to a map of observed erosion and deposi- 
tion for the dune. By comparing maps of the sediment 
budget (erosion and deposition) , it was possible to deter- 
mine qualitatively the accuracy of the wind tunnel simu- 
lation and the validity of the methods used in the experi- 
mental calculation of the sediment budget. 


peeceding page blank not ¥UMW 
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velocity of the wind increases with height, imparting hori- 
zontal momentum to the grains, so they travel a greater 
horizontal distance on the descent than on the ascent. 

The smallest grains may sometimes be expelled by 
larger grains and remain aloft for longer periods of time 
than the larger grains; those grains are said to be in 
suspension. Only a fraction of a percent of the total 
bulk of sand grains are transported in suspension [7] . 

The remaining 20-25% of particles are transported by sur- 
face creep, which is characterized by a rolling motion of 
grains on the surface. The largest grains generally move 
in surface creep. 

From reference 6, White et al. list several kinds of 
forces acting on an individual grain in saltation in a 
fluid medium: 

•V 

(a) the weight of the grain 

(b) drag, a force, tangential to the direction, of 
grain motion 

(c) lift, a force, normal to the direction of grain 
motion, caused by the pressure distribution on the particle’s 
surface 

(d) cohesional forces acting between individual 
grains 

(e) the force of overturning moment in shear flow 

(f) the Bassett force, which accounts for the effect 
of deviation in flow pattern from steady state 
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(g) the apparent mass force;, which accelerates the 
virtual mass of the particle relative to the ambient fluid 

(h) the force resulting from the effect of tempera- 
ture gradients in the flow field. 

For a more detailed description of these forces , see refer- 
ence 9. An additional force, resulting from the fact that 
grains in saltation are not isolated, but interact with 
grains on the surface, must also be taken into account; 

P. Nemenyi [8] calls this force the surface force. The 
surface force is any normal or tangential force exerted 
on a grain by any other grain, excluding cohesional forces. 

According to Hinze [9] , the Basset force and apparent 
mass force are only important if the density of the fluid 
is equal to or greater than the density of the particle. 
This is clearly not applicable to the case of sand in 
saltation, where the density of quartz grains is at least 
three orders of magnitude greater than the density of 
air [7]. From reference 6, the force of overturning mo- 
ment in a shear flow is of importance in calculating the 
minimum wind speed necessary to initiate saltation (fric- 
tion threshold speed), but is not important in calculating 
particle trajectory after initiation of motion. The ef- 
fect of the overturning moment force is discussed in 
greater detail in section 2.2. From reference 7, cohe- 
sional forces are not generally felt to be important for 
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sar,d grains in saltation, since most of the particles 
in saltation are too large to allow cohesional effects 
to be dominant. The temperature is assumed constant for 
the turbulent boundary layer, so forces resulting from 
the effect of temperature gradients will not be considered. 
The remaining forces can be divided into two classifica- 
tions; those responsible chiefly for initiation of the 
particle trajectory and those that affect the particle 
once it is in flight. 

For grains already in flight, Bagnold singles out 
the forces of gravity and fluid resistance acting in a 
direction opposite to the relative motion of the sand 
particle through the fluid as being the predominant forces 
acting on the grain [7] . 

Owen [5] gives the nonlinear equations of motion of 
the particles (neglecting lift) , once saltation has been 
initiated, as: 

my + mg + R(s)y/s =0 2.0.1 

mx - R(s) (U - x)/s = 0 2.0.2 

where the oncoming wind blows in the x-direction and y is 

. 21 ^ 

positive upwards. R is the drag and s — {y^ + (U - x) } 

From reference 6, the drag force, R, for a spherical parti- 
cle may be expressed as a function of the relative veloci- 
ty of the sand particle and the fluid and the coefficient 
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of drag, Cj^, in the following manner: 


R 


Cr^P 

D^a rel 



2.0.3 


where, for Reynold’s numbers less, than 100, 

=D = M W 

for 10 0 < Re < 2000, = iqO- 305865 (2 - logioRe) , 

and for Re > 2000, Cj^ = 0.4. Equations 2.0.1 and 2.0.2 

neglect lift as a force effective on particles after the 

initiation of saltation. It is unknown if lift is an 

important factor in determining the trajectories of sand 

grains; if an expression for lift, F , were included the 

lift force (to be defined later in this section) , would 

be subtracted from the left hand side of equation 2.0.1. 

The wind strength necessary to initiate saltation, 

u^ , is a function of the grain size (Cf. Fig. 2), 
t ' 

That a stronger threshold wind velocity is necessary 

to move larger grains is intuitively evident, but the 

fact that, for grains below a certain diameter, stronger 

winds are again necessary to initiate movement, requires 

some explanation. Below Dp Z 0.1 mm, the particles are 

too small to protrude above larger particles on the bed, 

and thus require a higher value of u^, to dislodge them 

t 

[ 10 ] . 
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Figure 2: Variation of the Threshold Velocity with Grain Size 

and Specific Gravity 

(Chepil, W. S. , "Dynamics of Wind Erosion II", Soil 
Science , 60: p.405, 1945) 



12 


There is some controversy concerning which forces 
cause grains to initiate saltation. Chepil, from reference 
2, discussed the forces initiating grain movement. He 
believed that, once saltation had been initiated, the 
impacts of saltating grains were responsible for continued 
saltation. He experimented with single sand grains on a 
smooth surface and noticed that a sand particle, after 
rolling a short distance, jumped vertically off the smooth 
surface, even though there were no obstructions. Combined 
Magnus and lift forces were suggested as a possible cause 
for the vertical rise. 

Prom reference 3, Chepil and Woodruff modified 
Chepil's earlier views on the importance of lift in the 
initiation of saltation. They said that the momentiom 
imparted to the grain by the lift force decreases rapid- 
ly with height and becomes hardly detectable a few grain 
diameter heights above the ground, considerably lower 
than the height to which many grains rise in saltation. 

Thus an additional vertical velocity, due to surface im- 
pact momentum transferred, must be imparted to the grain. 
From pressure measurements made on suspended spheres, 

Chepil and Woodruff conclude that lift cannot possibly 
be the sole factor involved in the vertical rise. The 
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angle at which a descending grain strikes a surface ob- 
struction would cause the rebound to be nearly vertical 
even if lift did not exist. (Similarly, the angle at 
which a grain in surface creep strikes another grain, 
coupled with the drag force, should cause the initial 
impulse to be vertical. 

Owen, from reference 5, believes that the vertical 
velocity of the particle is attributable to the geometry 
of grain impacts. He uses an order of magnitude analysis 
to find the ratios of lift and drag to the particle's 
weight and concludes that, for saltation in air, lift 
is at least an order of magnitude smaller than drag 
(assuming a spherical grain shape) . For saltation in 
water, where p/a is comparable to unity, lift plays an 
important role. (The ratio of the density of water, p, 
to a, the density of silicon dioxide, is 1:2.65, if 
buoyance effects are neglected) . 

White et al., from reference 6, assume that the sand 
particle velocity in the initial stages of lift-off is 
very small; they use that assumption as a justification 
for using Saffman's equations for lift on a spherical 
particle [11] . They maintain that, since saltation occurs 
in laminar flow [12] where there are no velocity compo- 
nents normal to the surface, the occurrence of saltation 
must be a result of a lift force. The laminar flow stud- 
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ied in the case above, however, occurred in a water stream 
where the density ratio between the sand particle and the 
fluid was comparable to unity. So although lift forces 
may be appreciable in the case of saltation in water, 
they cannot be assumed to be important to saltation in air 
The use of Saffman's equations to determine the lift force 
on a particle in saltation in air is unjustified, since 
one of the limiting conditions for using the equations 
is that the Reynold's number, Re = U^^j^^Dp/v , must be 
small compared to unity. If u* = 0.4 m/s is substituted 
for then, for a grain with Dp = 3 x 10"^* m, 

(v = 1.4 X 10“^ m^/s) , the Reynold's n\imber is 

8 . 6 . 

A numerical study was performed to calculate grain 
particle trajectories (neglecting lift) for different 
values of variables Dp, u* , z and V , the initial 
vertical velocity of the grain. Equations of motion, 

2.0.1 and 2.0.2, were numerically integrated, and values 
for X, y, X and y were calculated for different values 
of time, t. was made to vary as a function of the 

Reynold's number [13]. The height of rise and path 
length of the grains were tabulated for different initial 
conditions. Taking a typical example, it was found that 
a grain for which Dp - 0.3 mm in a flow where u^ = 0,4 m/s 
must have an initial vertical velocity of at least 2.5 m/s 
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in order to attain the height and downstream distance 
travelled found for the particle experimentally [14]. 

Thus Re for the particle is approximately equal to 50, 
which is obviously much greater than unity. 

Using the results from the computer program mentioned 
above, equation 2.0.3 was evaluated to find an approximate 
value for the drag on a particle within a few grain diam- 
eters* height of the ground. It was found that the drag 
on such a particle was about 5 x 10”^ Newtons. The weight 
of the particle, mg, is equal to 3.7 x 10”^ Newtons. The 
lift due to circulation around a spinning particle is 
given by: 

= p Ur 2.0.4 

L a 

where is the density of air, U is the velocity in the 
x-direction at a particular height, and r is the circula- 
tion. The lift, drag and weight were calculated for 
u* equal to 0 . 4 m/s, a particle diameter of 0.3 mm and 
equal to 0.08 ram. The circulation, r, is defined as 
follows: 

r = /U • ds 2.0,5 

where s is a line element on the spherical sand particle. 
The expression for the circulation must be integrated over 
the sphere, resulting in an equation for the lift force 
as follows : 

F — 2ttp UtoR cos^ ( — ) db 2.0.6 

L a m •' -R 2R„ 
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where R r the radius of the sphere is 0.15 mm, and oi is 
in 

the angular velocity of the particle. Integrating: 

2.0.7 

Values for U and uj are calculated using a logarithmic 

velocity profile (equation 2,1.2) and the fact that 

U = (dR . When these values , and values for p and R 
m m 

are inserted, the value for the lift force is found to 
be 6.84 X 10“® Newtons. Thus lift on a spherical particle 
is an order of magnitude less than either the drag or the 
weight. All studies that have been made of the forces on 
sand grains assume a spherical grain. Since, in 
particular, lift forces would vary greatly depending on 
the shape of the particle in saltation, an investigation 
ought to be made of the forces on different shaped sand 
grains . 

In conclusion, there are two forces that predominate 
in the saltation of a spherical sand grain once the grain 
has been ejected; the fluid resistance (drag) and gravi- 
ty. The initial vertical velocity of (spherical) parti- 
cle; ejection is provided chiefly by the momentum trans- 
fer that occurs during surface impacts between grains. 

The descending grain in saltation possesses momentum ac- 
quired from drag and previous impact. If it can be assumed 
that spherically-shaped sand grains are capable of initi- 
ating saltation, they must acquire an initial vertical 
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momentum from some source other than lift; one hypothesis 
is that they roll a short distance in surface creep as a 
result of drag, then rebound against surface obstructions 
vertically, beginning saltation. The nature and magnitude 
of lift forces on actual sand grains is unknown, but it 
is probable that lift plays an important role in the 
saltation of irregularly shaped grains and may not prove 
to be negligible, as it is when a spherical grain shape 
is assumed [15] , 
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2,1 Wind Velocity Profiles above a Sand Surface; 
Determination of Local Friction Velocities 
One of the first steps in the determination 
of erosion and deposition patterns at various locations 
on a barchan dune is finding the local friction veloci- 
ties at these locations, since the rate of sand flow is 
a function of the local friction velocity. Values 
for friction velocity may be determined from velocity 
profiles measured above the locations in question. 

Velocity profiles above a flat sandy surface have been 
measured by many experimentalists, and two types have 
been observed: profiles above stationary sand surfaces 

and profiles above drifting sand surfaces. 

From reference 7, Bagno Id writes the logarithmic 
velocity profile equation for a stationary sand surface 
(where saltation has not occurred) : 

U = 5.75 u*logio ~ 2.1.1 

'^o 

where U is the fluid velocity at height z above the sur- 
face with a roughness constant and friction velocity 
vA o/p^ , where is the surface shear stress and 
is the density of air. An alternate form of equation 
2 . 1. 1 is 
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where k, the von Karman constant is equal to 0.4. 

When the friction velocity is sufficiently large, 
the sand particles start to move. Bagnold defines 
two different threshold friction velocities initi- 
ating grain movement. The fluid thresho ,d is the 
friction velocity that the fluid must attain before 
the grains start to move without any other upwind inter- 
ference. The fluid threshold varies with the grain 
diameter as follows; 



where A is a constant coefficient which is 0.1 for salta- 
tion in air; Zingg [16] and Chepil [17] found the same 
value for A experimentally. a and p are the grain and 
air densities, respectively, and Dp is the particle diam- 
eter. If oncoming grains from upwind provide an initial 
source of kinetic energy, a lower threshold friction 
velocity is required to initiate saltation; this is 
the impact threshold. For grains where Dp > 0.25 mm 
(0.25 mm is the size of fine dune sand grains) , the impact 
is given by equation 2.1.3, where A = 0.08. Chepil, 
from reference 10, questions the relevance of adding kinetic 
energy to determine a separate threshold velocity. Recog- 


nizing the difficulty in determining visually the actual 
velocity at which saltation begins, he defines "minimal" 
and "maximal" threshold friction velocities. 

Once saltation is initiated, the velocity still varies 

s, 

as the log of the height, but the velocity rays pass 

through a new fixed focus at height k*, which is of the 

order of magnitude of the surface ripples formed by the 

sand flow [Cf . Figure 3] . The resulting ‘equation for 

the wind velocity U , where U is the wind velocity at a 

s s 

height z for the case of saltating sand, is 

where is the impact threshold velocity defined by 


u 




*t , k 
— In — 
< z. 


2.1.5 


is the impact threshold friction velocity defined 
in equation 2.1.3) . Experiments performed by Zingg indi- 
cate that, for wind over drifting sand surfaces, k is more 
closely approximated by ic = 0.375, but Horikawa and Shen 
[18] concluded that there is not sufficient evidence to 
deviate from a value of ic = 0 . 4 . 

For average dune sand of mixed grain size {Dp vary- 
ing from 0.1 mm to 1 mm), Bagnold estimates values of 
= 4 m/s and k' - 1 cm. These values are supported 
by Zingg' s experimentally determined approximations of 
k' = 10 Dp and = 20 Dp; k' and Dp are in millimeters, 
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and is given in miles per hour. Local values of the 

friction velocity are determined from the velocity pro- 
files at different locations on and around a barchan dune 
using equation 2.1.4 and substituting in Bagnold’s values 
for and k' . Values for U were obtained in the wind 
tunnel at a height z = 2.54 mm (for each location), which, 
when multiplied by a factor of 315 (used to scale field 
dimensions to wind tunnel dimensions) , is equal to a 
height of 0.8 m in the desert. The velocity is multiplied 
by a scaling factor, C, in order to normalize tunnel 
velocities to different assumed wind strengths at a height 
z = 0.8 m upstream with respect to the dune in the desert. 
This may be done because results from Chapter 4 show that 
the flow is Reynold's number independent for the values 
of Reynold's numbers at which the measurements were per- 
formed. From velocity profiles measured in the tunnel 
at a tunnel reference speed U_ = 10 m/s , the mean veloci- 
ty found upstream of the dune model at z =2.54 mm is 
4.5 m/s. Thus, in order to normalize all velocity values 
to an upstream velocity at that height of, for example, 

10 m/s, a scaling factor C = 2.22 would be used. 

Making the appropriate substitutions, the equation 
for the local friction velocity is 
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The method described above for determining the friction 
velocity assiimes a logarithmic profile between the lowest 
point of measurement (2.54 mm) and the point (U^,ln k') 
through which the velocity curves pass during saltation. 
For velocity profiles over a curved surface, the assump- 
tion is not strictly true and must be taken as only an 
approximation. 
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2.2 Parameters Affecting the Sediment Budget for the 

Barchan Dune 

From equations in section 2.1, an expression for the 
local friction velocity was derived from velocity profiles 
measured in the wind tunnel, where values were assigned to 
to k’ and U^. The local friction velocities thus deter- 
mined can now be used to model the sediment budget for 
the dune. For a flat sandy surface, the local rate of 
sand transfer, q, is a function of the local friction 
velocity and the grain diameter. 

From reference 7, Bagnold's equation for q, the mass 
of sand transported per unit time per unit width, as a 
function of u*', the friction velocity for moving sand 
over a flat surface, is: 



where c is a dimensionless, experimentally determined constant 
(C = 1.8 for naturally graded sand of the type found on desert 
dunes, and C = 1.5 for uniform sand grains where Dp = 0.25 mm). 
The particle diameter is denoted by Dp, D is a standard particle 
diameter (D = 0.25 mm) , and p is the density of air. 

Equation 2.2.1 is used to determine the rate of sand 
flow at the crest of the barchan model. The Crestline of 
the dune is the line at which the slope changes abruptly. 
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The slip face is in the lee of the crest and the leading 
edge is just upstream of it. In general, the flow separates 
at the crest and forms a reverse circulation at the slip 
face [19] . 

For the dune to be self-preserving, it must be in 
equilibrium; that is, all parts of the dune must have the 
same rate of translation in the direction of the mean wind 
and the height of the dune must remain constant. From 
reference 1, it can be shown that at the crest, the 
direction of sand transport closely follows the local 
wind [Cf . Appendix A]. In discussing the rate of sand 
flow over the dune, the emphasis up to this point has been 
placed on the effects of the local friction velocity on 
q. The barchan dune has a sloping surface; the effects 
of the slope angle, 9 (cf Figure 4) , on q must be investi- 
gated. Because of the surface slope, gravity no longer 
acts in a direction normal to the surface; particles 
initially moving will possess a downslope component of 
motion. Coarse particles travelling in surface creep are 
probably affected more than any other particles by surface 
slope. From reference 1, some particles very close to 
the dune base may be deflected by as much as 30® down- 
slope with respect to the direction of the surface wind. 
According to Bagnold (reference 7) , the coarsest particles 
generally determine ripple formation; thus surface slope 
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effects may be responsible for deviations measured between 
the local surface wind and the normal to the ripple strike 
at various locations. 


Bagno Id [20] modifies the equation for the role of 
sand transport as a function of friction velocity to take 
into account local variation in the gravity slope as fol- 
lows : 



P. 

g 



2 . 2.2 




where all variables are the same as in equation 2.2.1, a 
is the angle of repose of the grain on the bed, and 6 is 
the surface slope angle [Cf . Figure 4] . The slope angle 
also affects the threshold velocity, U^, of the sand 
particles. From reference 7, saltation occurs when 
F^/Fg > tana where F^ and F^ are the fluid and gravita- 
tional forces, respectively, on the grain. 

F^ = 0'pu*^Dp^ 2.2.3 

where B ' is an experimentally determined dimensionless co- 
efficient which depends on the turbulent fluctuations in 
the flow, the drag per unit area of the fluid on the grain 
and the height at which the drag force acts . 

Fg = J (a - p ) g -j 2.2.4 

and 
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At the initiation of particle movement. 


= tana = 


6 ' p u* ^ 

^ Dp(a - p)g 


2.2.5 


From reference 1, if surface slope effects are taken into 
account/ tana must be replaced by the function below; 


^ = /tan^a 


cos^e - sin^6 sin^Y 


COSY sine 

2 . 2.6 


where a and 6 are defined above and y is the angle between 
the surface wind and the surface gradient [Cf. Figure 5]. 

Equating 2.2.5 and 2,2.6, an expression for u* at the 
initiation of grain movement is easily derived; 

u^ = ' p ^ ^ [^/tan^a cos2^0 - sin^0 sin^y “ cosysine^ 

2.2.7 

From equation 2.1.2, an expression for the local threshold 
velocity is obtained and is inserted into equation 2.1.6 

to yield a new value of u* for each location. The adjusted 
u* values are used in equation 2.2.2 to obtain the local 
rate of sand transport. 

The orientation of the dune with respect to the mean 
wind affects both direction and quantity of sand transport 
at various locations on the dune. The equation used in 
Appendix A to find q for equilibrium at the crest includes 
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the angle, C, between the mean wind and Crestline. Values 
for £ are measured in the field. 

Two other factors that influence q are the lag effect 
and the effect of variation of grain diameter. When the 
local friction velocity at a given location, say, 
changes, the cimount of sand ejected into saltation changes 
also. Thus the amount of sand transported to a second 
location, say %2.i downstream of depends, to an extent, 
on wha' i akes place at £i. This "lag" effect should take 
place over a distance of the order of one path length of 
the grains in saltation, where the path length is the 
downstream distance travelled by a grain in a single jump. 
The path length varies from a few centimeters to a meter. 
Thus the concentration, q, of sand initially arriving at 
a particular location should be determined by what is 
happening at various distances upstream. The grain size, 
Dp, varies with different locations on the dune. As yet, 
this factor has not been taken into account in calculating 
the sediment budget on the dune; instead, a uniform grain 
size of 0.25 mm has been assumed. From equation 2 . 2 . 1 , 
however, it is clear that grain size affects the sand 
transport rate. To account for grain size effects, an 
accurate map of the distribution of Dp on a barchan dune 
would be required. 
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In a procedure more fully described in Chapter 4 , 
streamline maps of barchan #1 are drawn and "cells" of 
approximately equal area are marked off between stream- 
lines. Values for the parameters described above that 
affect q are found from field and wind tunnel measure- 
ments and are tabulated for locations at the center of 
each cell. A computer program provides adjusted values 
of q for each cell when different combinations of param- 
eters are taken into account. The effects of lag on q 
are computed for different assumed particle path lengths. 
The resulting erosion ‘md deposition rates are compared 
to erosion rates calculated for self-preservation of the 


sand dune. 



Chapter 3 


DESCRIPTION OF EQUIPMENT AND EXPERIMENTAL PROCEDURE 

3.0 Wind Tunnel Simulation of the Atmospheric Boundary 

Layer 

Techniques that have been developed to simulate the 
atmospheric boundary layer were adapted for the experiment 
to produce flow' over a desert surface [21] . The facility 
used for the simulation of a neutral desert boundary layer 
is an open return, subsonic wind tunnel with a test sec- 
tion 7.92 m long, 0.6 m wide and 0.61 m high (Cf. Figure 6). 
The air in the tunnel is initially forced through a con- 
traction with a ratio of 22:1. 

For a neutrally stable atmosphere, the desired mean 
velocity profile near the ground (where Coriolis effects 
may be neglected) may be written as 


U{z) 



3.0.1 


/ T 


where u* is the friction velocity ( U =y is the shear 

stress at the surface and p is the density of air) , ic is 
the von Karman constant ( =0.4) and is the roughness 
height [6,22]. A variable rod grid (rod diameter, d, is 
19 mm) produces, in conjunction with the appropriate rough- 
ness element fastened to the floor, the desired velocity 
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profile corresponding to the found for saltating sand 
[1] . The roughness element used is 20-3^5 grid silicon 
carbide sandpaper. From reference 21, two grids inserted 
at intervals upstream with respect to the rod grid pro- 
duce turbulence intensities similar to those found in the 
earth's atmosphere. Figure 7 is a plot of the turbulence 
intensity profiles. 

The scaling factor used in designing the experiment 
is 315:1. Thus the lowest 190 m of the atmospheric surface 
layer is represented by the tunnel height of 0.6 m and a 
2 ^ of 25.2 mm in the atmosphere is simulated by a of 
approximately 0.08 mm in the tunnel. A semi-logarithmic 
plot (Cf. Figure 8) of the mean velocity profiles taken 
at different intervals (x) downstream with respect to 
the rod grid shows good agreement with equation 3.0.1. 

The wind tunnel simulation of the atmospheric boundary 
layer compares well with naturally grown boundary layers 
in other wind tunnels [22,23,24,25]. 



Dimensionless Turbulence Intensity 
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3.1 The Sand Dune Models 

Velocity profiles are measured around three different 
barchan dune models, which are scale models of real dunes 
existing near the Salton Sea in California. In the field, 
the dunes were surveyed and contour maps were drawn [1] . 
Velocity profiles were measured, using cup anemometers, at 
various points on the dunes, with special emphasis placed 
along the crestlines. Also, maps were made showing the 
ripple marks, slope gradient, Crestline strike and slope 
curvature for each dune. 

The models , made from the contour maps , which will be 
labelled barchan #1, #2 and #3 for convenience, differ in 
some respects. Barchan #1 [Cf. Figure 9] is about 360 mm 
long and 310 mm wide (across the wing-tips) , barchan #2 
[Cf. Figure 10] is 430 mm long and 330 mm wide and barchan 
#3 [Cf. Figure 11] is 340 mm long and 200 mm wide. Barchans 
#1 and #2 possess a similar shape with symmetrical wings, 
indicating laterally uniform upwind sand transport, but 
barchan #3 i-s skewed (having one wing substantially dwarfed) , 
the result of an uneven sand supply upwind (from references 
6, 7). 

In the field, most of the measurements were made 
around the dune corresponding to barchan #1, so the most 
detailed set of profiles is determined for barchan #1 in 
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the tunnel. Barchan #1 also has symmetrical wings, and 
its dimensions are such that any effects from the warped 
walls of the tunnel did not interfere with velocity pro- 
files measured on the wings (the effects of the wind tun- 
nel walls were determined with a pitot tube; velocity 
readings were determined laterally at the position where 
the dune models were later inserted). Barchan #2 cannot 
be rotated through large angles (to measure effects of 
multidirectional winds) because it is so large that the 
tunnel's wall effects, which extend about 70 mm from 
each wall, do interfere with mean velocity readings on 
the wings . 

The dune models are carved from layers of balsa wood, 
with each layer corresponding to a height contour on the 
map. They are sanded to a smooth, continuous surface and 
spray-painted. The models are embedded into the floor 
roughness at a location 151 d (d = rod diameter) downstream 
of the rod grid. Barchan #1 is mounted on a disk that can 
be rotated through designated angles to simulate change of 
wind direction. Figure 12 shows the locations at which 
vertical velocity profiles are determinsrl on barchan #1. 

Line 7 is the center line of the dune and, in the case of 
winds impinging from 0°, is aligned with the axial direction 
of the tunnel and thus v;ith the freestream flow. 
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the tunnel. Barchan #1 also has symmetrical wings, and 
its dimensions are such that any effects from the 
walls of the tunnel did not interfere with velocity pro- 
files measured on the wings (the effects of the wind tun- 
nel walls were determined with a pitot tube; velocity 
readings were determined laterally at the position where 
the dune models were later inserted) . Barchan #2 cannot 
be rotated through large angles (to measure effects of 
multidirectional winds) because it is so large that the 
tunnel's wall effects, which extend about 70 mm from 
each wall, do interfere with mean velocity readings on 
the wings . 

The dune models are carved from layers of balsa wood, 
with each layer corresponding to a height contour on the 
map. They are sanded to a smooth, continuous surface and 
spray-painted. The models are embedded into the floor 
roughness at a location 151 d (d = rod diameter) downstream 
of the rod grid. Barchan #1 is mounted on a disk that can 
be rotated through designated angles to simulate change of 
wind direction. Figure 12 shows the locations at which 
vertical velocity profiles are determined on barchan #1. 
(The profiles were measured, for the most part, at the 
intersections of lines 1-13 with lines A, B and C.) The 
crest of the dune, which coincides v/ith line A in the area 
to the right of line 9, is indicated by the heavy dotted 
Point a occurs at the intersection of lines 1-13. 


line . 
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B, D', E*, F', G', H' and I' are points at which velocity 
profiles are measured in the lee of the model. Line 7 is 
the center line of the dune and, in the case of winds im- 
pinging from 0°, is aligned with the axial direction of 
the tunnel and thus with the freestream flow. 
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3.2 Velocity Measurement Apparatus and Procedure 

Mean velocity profiles and fluctuating components 
are measured with a miniature hotwire probe connected to 
a constant temperature anemometer. The overheat ratio 
for the hotwire is 1.6 times the probe resistance. The 
signal is fed through a linearizer and conditioned by 
an auxiliary unit with high and low pass filters and a 
nulling circuit which zeroes the mean voltage. D.C. vol- 
tage is measured using a digital voltmeter with a variable 
damping control which is kept at t = 10 seconds to average 
the u-velocities . 

The hotv/ire is positioned with a traversing system 
suspended from the tunnel ceiling, and it can be moved 
in the x, y and z directions; the traverse mechanism is 
accurate to 0.5 mm in all directions. Two kinds of minia- 
ture probes are used (Cf. Figure 13). The first type, 
which is designed with the sensitive element perpendicular 
to the prongs so that the probe support is horizontal and 
parallel to the x-axis , is used for measuring profiles on 
the dune's crest and on the lee side of the crest. The 
second type is a 45° probe (the sensitive element is fas- 
tened to the prongs at an angle of 45°) and is supported 
at an angle of 45° with respect to the z-direction so that 
the sensitive element is perpendicular to the wind; the 45° 






45® hotwire 





Figure 13: Perpendicular 


44 



45® Hotwires 
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probe allows accessibility to the upwind, slopes of the 
dune that would otherwise be inaccessible due to inter- 
ference from the probe support. 

A system was devised whereby hotwires could be cali- 
brated and the wind tunnel drift monitored without the 
insertion of a pitot tube. A manometer, having resolu- 
tion of 0.05 mm of water, is connected by rubber hoses 
to static pressure ports upwind of the contraction and 
at an x-location on the test section, close to the dune 
models. The manometer reads the differena© in static 
pressure, from which the mean wind velocity can be 
calculated using Bernoulli's eguation. With logarithmic 
and turbulence grids removed, the "reference" manometer 
reads the same values for wind speed as a pitot tube 
positioned at z^. With grids and roughness elements 
inserted, the relationship between the velocities determined 
by the pitot tube and the reference manometer is linear. 

For each profile measured, the tunnel drift was thus 
monitored and was restricted to 0. 4% of the mean velocity 
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3.3 The Use of Flow Visualization as a Wind Direction 

Indicator 

Hotwire measurements of profiles around dune models 
are sufficient to provide information about relative wind 
speeds at various spatial locations. However, only in 
one case is the hotwire used to directly indicate flow 
direction; this case will be discussed later in this 
section. In all other eases, the technique of flow visu- 
alization is employed to determine directional trends . 

Several methods of flow visualization were investi- 
gated; the method of tufting proved to be the most success- 
ful. Thin strips of tissue paper, the "tufts," are glued 
at one end to the models being studied, attached in such 
a way as to have no initial directional bias. When the 
wind blows on the models, the tufts align themselves with 
wind direction, thus indicating streamlines on the dunes. 
Photographs are taken of the models at different exposure 
times to determine both "instantaneous" and "averaged" 
wind trends. A sample of tuft photography for barchan #1 
is shown in. Figure 14; 14(a) shows an instantaneous tuft 
alignment with an exposure time of 1/2 second and 14 (b) , 
taken at an exposure of 2 seconds, illustrates the average 
wind flow. It is from photographs of the averaged cases § 

that streamline maps of the flow over the dune are drawn. 
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At points near the base of the dune on the windward 
side, the tufts normally remain in one position, not 
always indicative of flow direction, because the wind at 
those points is not strong enough to move the tufts effec- 
tively. Clearly, for these points another method for find- 
ing wind direction is needed. it is in such a case that 
the hotwire is used as a direction indicator. The hotwire 
has an angular cosine sensitivity to the flow and may be 
rotated to find flow direction [27]. A 45“ hotwire does 
not need to be rotated; it may be inserted into a stationary 
probe support in two different orientations: the first 

at an angle of 45® with respect to the x-axis in the posi- 
tive y direction; the second at an angle of 45“ in the 
negative y direction (Cf . Figure 15) . Using the equation: 

E2 — E]^ 

E2 4- El 

wheire E^ and E 2 are voltages (corresponding to velocity 
magnitudes) recorded at the two orientations, it is easy 
to compute the resultant angle, 0', at which the wind is 
blowing with respect to the x-axis. The angles thus calcu- 
lated contribute to the streamline maps of the model dunes. 
In the next chapter, the resulting streamline maps and 
their importance in computing dune migration wi ) 1 be ex- 
plored. 
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Chapter 4 


EXPERIMENTAL RESULTS 

4.0 Velocity Profiles around a Barchan Sand Dune 

In. order to determine local friction velocities 
which are needed for the analysis of the sediment budget 
of a barchan dune, it was found necessary to measure an 
extensive set of velocity prof iles . Hotwire measurements 
of U were made up to heights far greater than was necessary 
to determine u* values; one of the additional purposes of , 
the measurements is to determine flow patterns of the wind 
above the saltation layer of the dune. Profile measure- 
ments provide the answers to basic questions about the 
interaction between the mean wind and the barchan. 

From a study of these velocity profiles, it is possible 
to determine certain flow trends; for example: the height 

at which the flow resumes a logarithmic profile undisturbed 
by the dune is approximately equal to six time the maximum 
dune height. A region of reverse flow in the lee of the 
dune crest is detectable with a hotwire. Changing the 
speed of the mean velocity of the wind has no effect on 
the velocity profiles other than a scaling effect; if the 
mean wind speed is doubled, for example, each of the veloci- 
ties in the profiles is also doubled, to within 4% accuracy. 
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Figure 17(a); Velocity Profiles Along Line 7 Fo 
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Figure 16: Velocity Profiles in the Lee of Barchan #1 

(See Figure 12 for locations of positions 
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at a decreasing rate as the crest is approached, so that 
the profiles converge and become logarithmic at z = 160 mm. 
Profile a, in the lee of the dune, possesses a very low 
velocity for z < 10 mm which actually decreases initially 
as z = 10 mm is approached, then greatly increases to con- 
verge with the other profiles. The hotwire is not able to 
adequately interpret reverse flows, so it can only meaning- 
fully be used at a to determine the location of this "dead" 
zone in the lee of the crest, and not its characteristics. 

The dead zone is seen to be approximately 10 mm high and 
extends no more than about 90 mm downwind from a. 

Similarly, profiles measured on the wings of the dune 
model were studied to determine flow trends. For conve- 
nience, the wing on which lines 8-13 are drawn is called 
wing A and wing B is the wing on which lines 1-6 lie. On 
wing A, between height contours 6 , 3 mm and 9 . 5 mm [Cf. Fig- 
ure 18], the wind attains its greatest near-surface (C' = 

2.54 mm) velocity at line 11, which is exactly perpendicu- 
lar to the oncoming mean wind . Thus , the flow experiences 
an acceleration along that location of the wing and a 
deceleration as the sheltered portion of the wing (lines 12 
and 13) is approached. Along the Crestline of wing A, the 
velocity at C' =2.54 mm decreases downstream; this is to 
be expected, because the height z of the Crestline decreases 
downstream also. A relatively high velocity (at = 2.54mm) 
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Figure 17 (b) : Velocity Profiles Along Line 7 for 

Barchan #1 
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occurs at the wing tip, indicating a convergence of flow 
at that point [Cf . Figure 26] . For all cases except at the 
tip, the velocity increases from the lower contours of 
the dune to the Crestline. All of the profiles measured 
on wing A converge at z = 60 mm (scales to a field height 
z ~ 19 m, approximately three times the maximum dune height) 
with each other. In general the profiles at the crest con- 
verge with each other at z = 20 to z = 25 mm. On both wings, 
between contours 1.27 and 1.59, the local wind is deflected 
slightly away from the wing tips and towards the center of 
the crest. 

Analyzing profiles of wing B, one notices that be- 
tween contours 0.63 and 0.95, the maximum near-surface 
ve locity occurs at line 4 , which is not perpendicuJ.ar to 
the oncoming wind, rather than at line 3, which is. The 
probable explanation for this is that at line 3 the wing 
is indented slightly, which has the effect of protecting 
it from the wind; at line 2, the wind accelerates again, 
where the wing protrudes outward. The profiles on wing B 
converge at z = 0.1 m, with the crestal profiles being the 
last to converge with the lower ones; the crestal profiles 
converge with each other, however, at z = 30 mm. 

A set of velocity profiles was measurei., for a tunnel 
rerference speec 0^^ = 10 m/s in the lee of the dune at six 
locations laterally positioned at a distance 210 mm down- 
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Figure 18: Velocities at ' = 2.54 mm. Superimposed on Height 

Contours for Barchan #1 
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occurs at the wing tip, indicating a convergence of flow 
at that point [Cf . Figure 26] . For all cases except at the 
tip, the velocity increases from the lower contours of 
the dune to the Crestline. All of the profiles measured 
on wing A converge at z = 60 mm (scales to a field height 
z = 19 m, approximately three times the maximum dune height) 
with each other. In general, the profiles at the crest con- 
verge with each other at z=20toz=25 mm. On both wings, 
between contours 12.7 and 15.9, the local wind is deflected 
slightly away from the wing tips and towards the center of 
the crest. 

Analyzing the profiles of wing B, one notices that be- 
tween contours 6.3 and 9.5 the maximum near-surface 
velocity occurs at line 4 , which is not perpendicular to 
the oncoming wind, rather than at line 3, which is. The 
probable explanation for this is that at line 3 the wing 
is indented slightly, which has the effect of protecting 
it from the wind; at line 2, the wind accelerates again, 
where the wing protrudes outward. The profiles on wing B 
converge at z = 0.1 m, with the crestal profiles being the 
last to converge with the lower ones; the crestal profiles 
converge with each other, however, at z = 30 mm. 

A set of velocity profiles was measured, for a tunnel 
reference speec = 10 m/s in the lee of the dune at six 

XV 

locations laterally positioned at a distance 210 mm down- 
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stream with respect to a [Cf. Figure 16] . Profiles E', 

F', G' and H' converge and begin to become logp,:ithmic 
at z - 95 mm; profiles D' and I* do not converge until 
a height z = 210 mm is attained. Velocities recorded 
for the profiles at a height z = 2.54 mm are studied in 

order to determine a flow pattern for the dune. 

\ 

A rough pattern of "streamlines" is drawn from the 
measurements [Cf . Figure 19] ; the pattern resulting at 
z = 2.54 mm coincides closely with Allen's streamlines 
in the lee of a barchan dune [19] . 

The effect on velocity profiles of changing U , the 
tunnel reference speed, was investigated. Velocity pro- 
files were measured at key points on the dune (7A, 13 
crest and 2 crest) and at a point 10 cm upwind of the 
dune at tunnel speeds of 5, 6, 7, 8 and 10 m/s. A scal- 
ing factor by which velocity profiles could be normalized 
to a single reference speed (10 m/s) was sought. The 
value of predictability of profiles for different on- 
coming wind speeds is that it can be used to cut down 
considerably on the number of measurements necessary to 
determine flow around the barchan dune. 

The most obvious method of determining scaling factors 
was investigated first and proved to be fruitful. Each 
point in the velocity profile is multipled by a scaling 
factor equal to U„/U where U_ = 10 m/s and U. is the 
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tunnel reference speed for the profile. The sctled pro- 
file is plotted on a graph with a profile measured at 
the same location for U = 10 m/s. The degree of collapse 
may thus be evaluated, or the velocities may be individual- 
ly compared to find a percentage difference between them. 
Figure 20 shows sample velocity profiles measured at 
7A at tunnel speeds of 5 and 10 m/s; the profiles are 
nondimensionalized by dividing by the appropriate reference 
speed. The two profiles coincide closely with an average 
error of 1.3%; most of the difference occurs close to 
the dune's surface. The other profiles yielded similar 
results with an overall average percentage difference of 
4%. Later in this section, the source and magnitude of 

the errors wilx be presented- It will be shown that 
the average percentage error in determining velocity 
profiles is about 3.8%, so an error of 4% may 
be said to lie within experimental error. In conclusion, 
the velocity profiles may be scaled in the manner described 
above to obtain a representation of the profiles at some 
desired wind velocity. Reynold's number effects are 
thus found to be negligible for wind tunnel speeds used 
in the experiment. 

An experiment was performed to determine the 
extent to which the logarithmic boundary layer artificially 
produced in the tunnel affects velocity profiles around the 
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dune model.. Profiles were measured at certain key loca- 
tions Oil the dune with the turbulence and logar5 thmic 
grids and the roughness element removed. Removing the 
grids and roughness element resulted in the formation of a 
naturally grown boundary layer in the wind tunnel very close 
to the tunnel floor; for convenience, this flow shall be 
called "uniform," since the flow above the boundary layer 
is uniform, up to the heights at which measurements were 
made. The profiles thus determined were compared to profiles 
measured at the same positions in the logarithmic boundary 
layer. Figure 21 is a sample comparison. It is evident 
that the two profiles are dissimilar. At the point closest 
to the surface of the dune, where the dune affects the flow 
to the greatest extent, the percentage difference between 
the two velocities is 5.7%. The difference increases to 
a maximum of 12.11% at z/z^ = 1000, then drops as the 
profiles converge at z/z^ = 2750, where the logarithmic 
velocity at that height becomes equal to the free-stream 
velocity. The two profiles separate again as z increases 
beyond z/z^ = 2750.. From the results above and similar 
results obtained from other profiles, it is established 
that a simulated boundary layer is necessary if an ac- 
curate representation of wind flow over a barchan dune 
in the atmosphere is to be obtained, provided that the 
atmospheric boundary layer is really logarithmic. 



Nondimensionaiized Mean Velocity ^ 



z 


^0 

Nondimensionaiized Height Above Tunnel Floor 

Figure 21: Comparison Between Free-Stream and Logarithmic 

Velocity Profiles 
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A brief investigation was made into the effects on 
velocity profiles of changing the angle of the v.-ncoming 
wind with respect to the dune. The dune was rotated 
(horizontally, about the axis at point a), through angles 
of -30° to 30° with respect to the mean wind, at 10° inter- 
vals. Profiles were determined, for U = 10 m/s , at posi- 

K 

tions 13 crest, 7A and 2 crest. Figures 22, 23, 24 and 25 
show the resulting velocity profiles for positions 2 crest 
and 13 crest. 

At an angle of 30°, wing A and, thus, 13 crest, is 
exposed more directly to the oncoming wind and wing B is 
partially sheltered. The "sheltering" effect can be ob- 
served in Figures 22-25. Velocities near the surface of 
13 crest are highest when the mean wind impinges at 30° 
and become progressively lower as the wind is rotated to 
an angle of -30°. The opposite effect, as expected, occurs 
for wing B. The velocity profiles converge at = 3000. 

The effect of changing the angle of attack of the wind on 
the velocity profiles and on the sediment budget of the 
barchan is one that demands more research. 

The predominant source of error in the measurement of 
mean velocity profiles above the model occur in the posi- 
tioning of the hotwire at locations close to the dune's 
surface. The distance between the sensitive element of 
the probe and the desired location on the model was 
measured with a vernier accurate to 0.01 cm, but some 
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difficulty was experienced, due to the extreme delicacy 
of the probe, in getting the vernier close enc-<gn to 
the sensitive element to make a reliable reading. Also, 
under the influence of the oncoming wind, the sensitive 
element often moved slightly from its initial position, 
requiring an additional adjustment of the hotwire’ s 
height, which had. to be estimated visually. Assioming a 
displacement of 1 mm for the hotwire from the desired 
vertical location at each position, the average percentage 
of error for a complete velocity profile is 2.2%. (This 
error was calculated using equation 2.1.2 and varying 
z by 1 mm) . The error from hotwire displacement varies 
from 15%, at C' = 2.54 mm (where is the vertical dis- 
tance) of the sensitive element from the surface where the 
profile is determined), to 0.1%, at C' = 230 mm. 

At locations close to the surface, particularly in 
the area of the crest and wings of the model, the turbu- 
lence intensities were quite high. The damping circuit 
on the voltmeter used to determine the mean velocity was 
sometimes not adequate to deal with turbulent fluctuations, 
and an estimate had to be made of the voltage, which on 
occasion fluctuated by as much as 2% of the mean velocity. 
Tunnel drift and hotwire calibration together accounted 
for about 0.2% error in all velocity profiles. Totaling 
percentages of the sources of error listed above, and 
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assigning an average value of 1.4% to turbulent fluc- 
tuations, the average error resulting from velocity pro- 
file measurements is equal to 3,8% of the mean velocity 
in the wind tunnel. 
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4 . 1 Local Friction Velocities on the Barchan Dune 

The method used for determining values for the local 
friction velocity, u^, on the barchan dune is discussed 
in section 2.1. Equation 2.1.6 yields values of friction 
velocity as a function of the velocity at a height z = 

2.54 mm in the tunnel, equivalent to a height of 0.8 m 
in the desert. Using the C = 2.22 in equation 2.1.6, 
friction velocity values are calculated for each location 
at which velocity profiles are determined on barchan #1. 
Figure 26 is a contour map of the local friction velocities 
drawn from a map of calculated u^ values . 
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4.2 Streamline Mapping of Barchan #1 

In determining the sediment budget for baichan #1, 
it is necessary to construct accurate streamlines for the 
dune. The streamlines are subsequently divided into 
"cells" of equal length, and values for erosion or deposi- 
tion of sand are computed for the center of each cell. 

Three different methods are employed in determining 
streamlines for barchan #1. The first streamline map 
[Cf. Figure 27] was drawn from a combination of tuft photo- 
graphs and directional measurements made at the base of 
the model with a 45° hotwire proble, using the method out- 
lined in section 3.3. The tufts in the photographs are 
aligned with the local wind; streamlines may be interpolated 
between the tufts from points slightly upwind of the dune 
base to points on the crest. Local wind vectors deter- 
mined with the hotwire at the dune ' s base are added to 
the tuft streamlines to form a complete streamline map of 
the model. 

Measurements made in the field for barchan #1 provide 
the second method for determining local wind streamlines ^ 
for the dune [1]. Stakes were implanted on the dune's 
surface and the directions of the scour marks resulting 
from wind action were measured using a compass. The 
"vectors" determined by the scour marks were measured 
for two different mean wind orientations, impinging from 
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N73E and from S81E, The desired local wind vectors are 
for the case where the mean wind impinges from lv32E (the 
orientation of the dune model in the wind tunnel) , so a 
"weighted mean" value for local wind direction at a mean 
wind of N82E is calculated from values obtained in the 
field. Streamlines are drawn from the direction indicators 
thus derived in combination with vectors obtained in the 
wind tunnel with. the hotwire at the dune’s base. Figure 28 
is a streamline map determined using the second method for 
finding local wind vectors. 

The third method for determining "streamlines" for 
a barchan dune assr*'es a local wind always perpendicular 
to local ripple marks. From reference 2 ., ripple strike 
nseasurements were made in the field for barchan #1 at 
the two wind orientations described in the previous para- 
graph, and vectors for the local wind direction for an 
impinging mean wind at N82E are interpolated from the 
two measured set;'> of vectors. Figure 29 is a map drawn 
from the local wind vectors derived from the ripple strike 
measurements; this map will be called a "ripple-line" 
map, because the lines on it are not streamlines. The sedi 
ment budget for the dune is calculated for the ripple-lines 
in Figure 29 in the same manner as for Figure 28; the 
paths determined’ by the ripple-lines are divided into cells 
of equal length,; and erosion and deposition rates are de- 
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termined for each cell. In Figure 29, certain ripple- 
lines on the wings branch off from others to fotni two 
separate ripple-lines. This occurs where the ripple lines 
are dense, indicating a strong convergence of local wind. 

Patterns of erosion and deposition found in section 
4.3 for the barchan dune, using the streamline and ripple 
line maps in Figures 28 and 29, will provide a comparison 
between the two "streamlining" methods. However, com- 
parisons may also be made among all three maps by exam- 
ining the figures. A comparison of the tuft streamlines 
and streamlines determined in Figure 28 from field mea- 
surements reveals many similarities between the two. 

Since the hotwire method described j.n section 3.3 was 
used for both maps to find directional trends at the dune ' s 
base, the streamlines in Figures 27 and 28 are initially 
identical. Both sets of streamlines diverge initially on 
the main section of the frontal slope* upwind of the wings, 
then converge at the crest. Near the tip of wing A, both 
streamline maps indicate a strong convergence of flow 
toward the tip, rather than toward the crest,* on wing B, 
both maps show a flow trend towards the crest for the local 
wind, rather than towards the wing tip. 

The principle difference between the streamlines in 
Figures 27 and 28 occurs on the sections of the wings 
close to the main central portion of the dune. Figure 27 
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shows slightly greater convergence of local flow towards 
the wing tips than Figure 28 does (where "convc/gence , " 
in this case, implies a bottleneck effect) ; the local wind 
in Figure 28 shows the greater flow tendency towards the 
crest. Because the two maps are so similar, only one, 
Figure 28, was selected for sediment budget analysis in 
section 4.3 ' ^ 

Comparing Figures 27 and 28 to Figure 29, it is im- 
mediately evident that ripple-lines show greater conver- 
gence of flow towards the wings than the streamlines do. 
Streamlines following the measured local wind directions 
indicate a comparatively stronger flow towards the cen- 
tral crest area of the dune. 

Flow lines for the maps in Figures 27, 28 and 29 are 
drawn by hand from photographs or draW-ings containing only 
the short, vector-like direction indicators determined 
in the field and wind tunnel. Some error must result from 
mis judgment in the interpolation of lines between the di- 
rection indicators. In an effort to minimize human error, 
the streamline maps were drawn independently several times 
by two different people, until the fihal copies were deemed 
suitable by both. 
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4.3 Mapping the Sediment Budget for Barchan #1 

The final stage of experimental results is ^ compari^ 
son of contour maps for the sediment budget of barchan #1. 
Contours obtained using various combinations of the equa- 
tions outlined in section 2 . 2 and the results of the 
streamline and ripple-line maps in section 4.2 are compared 
with each other and with predicted and observed contours 
of erosion and deposition for the dune. All erosion and 
deposition contours are normalized with respect to the 
rate of forward motion of the dune. the erosion at 

a point 1 on the dune's surface, is equal to Ah/Ax where 
Ah is the decrease in the height of the surface at point 
1 when the dune has migrated a distance Ax. Thus E and 
D, the deposition, are nondimens ional ? E is assigned nega- 
tive values because h decreases during erosion, and D is 
assigned positive values since h increases with deposition. 

Figures 30 and 31 are erosion and deposition contour 
maps drawn using the streamlines of Figure 28. Values 
for u* used to obtain q, the sediment transport rate, are 
obtained from velocity profiles measured in the wind tun- 
nel. The two sets of contours differ only by the slope 
factor used to determine q. Figure 30 is drawn using 
Bagnold's equation, as expressed by equation 2.2.1, for 
sand transport, neglecting the slope effects of the dune. 
Figure 31 results from using Bagnold's equation, as ex- 
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Figure 30: Erosiorx and Deposition Contours for Barchan #1 

( Neglecting Slope Effects, Using Streamlines from 
Local Wind ) The values represent dimensionless 
erosion and deposition, where negative values indicate 
erosion and positive values indicate deposition. 



Figure 31: Erosion and Deposition Contours for Barchan #1 

( Including Slope Effects, Using Streamlines from 
Figure 12 ) , 
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pressed in equation 2.2.2. Contours for the tangent of the 
slope angle of the barchan, from reference 1, aie given in 
Figure 32. From a visual comparison of the two figures, 
many similarities are evident. The contours on wing A 
in both maps are nearly identical. The portion of wing 
B near the central section of the’ dune in Figure 31 ex- 
periences less erosion than the corresponding location 
in Figure 30 does. In Figure 31, the main section of 
the dune experiences less erosion and the central crest 
section more erosion than in Figure 30. In general, the 
two maps show erosion in the central section of the dune, 
and on the wing areas closest to the central section. 
Deposition on the dune (as opposed to deposition in the 
lee of the crest) occurs on the wing tips. 

Figure 33 is a contour map obtained from values of 
erosion and deposition found for the ripple-lines drawn 
from ripple marks in Figure 29. Again, Figure 33 is 
basically similar to Figures 30 and 31, where q is de- 
termined by the local wind. Erosion occurs on the main 
section of the dune and deposition takes place on the 
wing- tips and on a large section of the crest of wing A. 

In section 4.2, it was noted that the ripple-lines of 
Figure 29 indicate a strong convergence of flow along the 
wings; in Figure 33, deposition does occur on a greater 
portion of the wings than in the other two figures , but 
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Figure 32: 




Figure 33: Erosion and Deposition Contours for Barchan #1 

(Including Slope Effects) Using Ripple-lines 
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there is relatively much less deposition on the crest of 
wing B. The erosion contours in Figure 33 on che main 
section shift towards the crest and show a higher erosion 
rate than either of the other sets of contours do. 

Figure 34, from reference 1, shows sediment budget 
contours for the dune drawn from measurements made of 
stake exposure and burial over a two week period. The 
pattern of erosion on the dune's central section differs 
slightly from all the other cases discussed so far, in 
that the area where the greatest amount of erosion occurs 
is closer to the crest than in the other cases. There is 
more erosion on the tip of wing B and more deposition on 
the crest of wing A than in the other figures, but the 
erosion on wing A is similar in extent and location to 
that found in Figures 30 and 31. 

Figure 35, from reference 1, is a map of the predicted 
erosion and deposition contours for barchan #1, if dune 
equilibrium is preserved [Cf . Appendix A] . Comparing the 
contours in Figures 30, 31, 33 and 34 with the predicted 
contours, it is evident that the sediment budget calculated 
for streamlines following the local wind bears the greatest 
resemblance to the predicted sediment budget. Erosion rates 
on the central portion of the dune are similar to those 
found for Figures 31 and 34. On wing B, the deposition 
on the crest for the predicted case is greater than for 



Figure 34: Erosion and Deposition Contours for Barchan #1 

(From Field Observations over a Two-Week Period of 
Stake Exposure and Burial, Reference 1) 


Figure 35: Predicted Erosion and Deposition Rates for Barchan 

#1 (Reference 1) 
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most cases and very similar to the deposition occurring 
at that location in Figure 34; the erosion on wing B is 
similar to that indicated in Figures 30 and 31. The 
erosion occurring on wing A in Figure 35 also takes place 
in Figure 30; deposition on wing A closely resembles the 
pattern found for Figures 30 and 31. 

In this section, patterns for the sediment budget 
discussed have taken into account only a few of the cases 
that were investigated. From reference 1, the factors 
of grain size, delay time and initial saturation of the 
oncoming wind discussed in Chapter 2 were varied; none 
of these factors appeared to have a profound effect 
on the sediment budget. The pattern of erosion and de- 
position found for streamlines following the local wind 
closely resembles the predicted pattern on the wings, but 
differs somewhat from the predicted pattern on the main 
section of the dune. 

Contour maps of the sediment budget for the barchan 
dune are the cumulative result of maps determined for the 
streamlines, slope angle and friction velocities on the 
model. Because of this, they are subject to all of the 
sources of error listed in this chapter so far, as well 
as some new sources of error resulting from the drawing 
of the erosion and deposition contours themeselves. Con- 
tours for the friction velocity in Figure 10 (drawn by hand 
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and subject to hxoman error) , are rather widely spaced in 
some locations, particularly on the leading edge of the 
dune. Thus, the values for u* found for the different 
cells is subject to additional error. Streamlines drawn 
from field measurements may have been determined when the 
dune was not in a state of equilibrium, a fact which could 
be the reason for the discrepancy between the predicted 
and "measured" erosion patterns on the central section 
of the dune. 

Although the sediment budget for the dune in the case 
of streamlines following the local wind direction bears 
a close resemblance to the predicted sediment budget for 
dune equilibrium, there is also great similarity between 
the erosion patterns on the central dune section for the 
ripple mark and predicted cases. Further investigations 
are required to determine factors needed in predicting 
the sediment budget; a model combining streamlines obtained 
from the local wind and rippJ.e marks may produce erosion 
and deposition patterns more similar to the predicted ones. 
In conclusion, there exists a promising indication that 
it is possible to predict the rate and direction of sand 
flow around obstacles Using the techniques of wind tiinnel 
modelling of conditions in the field and hotwire measure- 
ment of velocity profiles. 



Chapter 5 


SUMMARY 

5.0 The Interaction of Unidirectional ^inds with a Barchan 

Sand Dune 

An experimental investigation was made of the flow 
patterns around a model barchan sand dune which was inserted 
into a wind tunnel simulation of the Earth's desert boun- 
dary layer. A listing follows of the principal results 
achieved. 

1. The desert boundary layer of the Earth was success- 
fully simulated, with a roughness value, z^, equal to 0.08 
mm, 

2. Velocity profiles measured with a hotwire anemom- 
eter provided insight about flow trends around a model 
barchan. 

3. A method for scaling velocity profiles to predict 
the effects of changing the speed of the mean oncoming wind 
was determined. 

4. A comparison of velocity profiles measured on the 
dune in uniform and logarithmic flows proved the necessity 
of utilizing a logarithmic boundary layer simulation in 
determining flow patterns around a barchan dune. 
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5. A brief investigation was made into the effects 
on velocity profiles of changing the angle of Lhe oncoming, 
mean wind, and resulting changes in the flow trends were 
noted . 

6. Streamline mapping using flow visualization methods 
produced results similar to streamline maps of the barchan 
dune obtained from field measurements. 

7. A theoretical investigation was made into the 
forces acting on a sand particle in saltation. 

8 . Values for the local friction velocity at various 

4 

points on the dune were derived from velocity profiles, 
using Bagnold's equation as stated in eq. 2.1.6. 

9 . Experimentally determined maps of the sediment 
budget for the barchan dune are compared to each other 
and to a map of the sediment budget determined theoreti- 
cally for dune equilibrium. 
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Appendix A 


COMPARISON OF THE RATE OF SAND TRANSPORT .^'iASURED 
AND CALCULATED TO SATISFY DUNE EQUILIBRIUM AT 
THE CREST OF BARCHAN #1 (FROM REFERENCE 1) 

The requirement for equilibrium of the dune ' s crest 
under the action of unidirectional wind is that the rate 
of translation of the crest be uniform at all points and 
that the crest height remain constant. The volume rate 
of delivery of sand per unit length, I, of the crest, is: 

q sin (?-«') 

where g is the angle between the mean wind and the Crest- 
line and a ' is the angle between the transport direction 
and the mean wind. 

The volume rate of building perpendicular to the 
crest is 

H & 

^ dt 

where y is the horizontal direction perpendicular to the 
Crestline and H is the height of the Crestline above the 
base of the dune. Because of curvature, the crest height 
must be replaced with an effective height, H^. At equi- 


librirun: 


He ^ sin? = q sin(? - a') 


or 
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He sing 
sinlc - ot * ) 


where k' = dx/clt., a constant on all parts of the Crestline 
Figures 36 and 37 show k* as a function of position 
on the Crestline/ for q measured in the field and in the 
wind tunnel (during the course of this study) along the 
crest, respectively. The value for k' shows some scatter 
due to experimental inaccuracies, but a good order of mag- 
nitude estimate of k' may be obtained from the graphs. 







wing A center wing 

Figure 37: k' vs. Position on Crestline for Barchan #1 (From Wind Tunnel Measure 

ments of q) 



